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Abstract 

Using e expansion proposed in [1] we calculate density correlation function of the degenerate 
Fermi gas at infinite scattering length to next-to-leading order in e for excitation energies below 
quasi particle threshold. An expression for the low energy dynamic structure function is derived. 
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I. INTRODUCTION 



Cold dilute gas made of two Fermion species (which will be labeled as |, I) with two body 
interactions characterized by 

na 3 > 1, nrl < 1, (1) 

where n is the number density, a is the s wave scattering length and ro is the effective range 
of the potential, named "unitary Fermi gas" is one of the simplest examples of a strongly 
interacting Fermi system. Yet, unitary Fermi gas is only a slight idealization of dilute 
neutron matter (a ~ — 18/m, ro — 2.6/m in the 1 So chanel) that may exist in the crust of 
a neutron star. Also, it has been studied extensively in the experiments on trapped cold 
atomic gases with tunable interactions using the Feshbach resonance technique [2-8] . From 
the theoretical point of view, unitary Fermi gas lacks any intrinsic scale and is expected to 
exhibit universal properties. At the same time theoretical description of unitary Fermi gas 
is difficult due to apparent absence of a small dimensionless parameter which could be used 
in a perturbative expansion. 

Inspired by ideas of Nussinov and Nussinov [9], an analytical technique similar to the e 
expansion in the theory of critical phenomena has been recently proposed [1]. The idea is 
to describe 3 dimensional unitary Fermi gas using perturbation theory in the dimensionality 
of space, d, around d — 4 and d = 2. The rational for such an approach is that in d = 2, 4 
unitary Fermi gas simplifies 

• to the non-interacting Fermi gas in d — 2, 

• to the Bose-Einstein condensate of non-interacting spin dimers in d = 4 

[9, 10]. Comparison with the results of Monte Carlo simulations suggests that already at 
the next-to-leading order (NLO) in e the expansion around d — 4 may be a useful tool in 
the description of this system [1]. While subsequent investigation raised doubts about the 
convergence of the series in e [11], one obvious way to check usefulness of the technique is 
to make predictions for various observable quantities at NLO. If the trend suggested by the 
initial NLO results of [1, 10, 12] were to hold, it would serve as a strong encouragement for 
further investigation. 

So, as a step in this direction, in this paper we report calculation of the density correlation 
function of the degenerate unitary Fermi gas in d = 4 — e spatial dimensions to the NLO in e 
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for energy below quasiparticle threshold. Expressions for the low energy dynamic structure 
function and static structure factor are quoted. The simplest way to make predictions for 
d = 3 is to set e = 1 in the NLO expressions (15,17,18). 



II. THE LAGRANGIAN AND DEFINITIONS 

Due to universality of unitary Fermi gas, we choose to work with local four-Fermi in- 
teraction with a coupling constant Co tuned to reproduce infinite scattering length. The 
Lagrangian density of the unitary Fermi gas in the presence of the source coupled to density, 
p(x), is (here and below H — 1) 



C = ifi ^id t + ^ + p + p(x)^j 4> + c n T n h 



(2) 



where ip = (ipi,ipi) T , is spin-1/2 fermion field, and rii = ipj ipi, with i = { f, j } are particle 
density operators, p is the particle number chemical potential which in this calculation is 
assumed to be the same for both pairing species. The theory is defined in d = 4 — e spatial 
dimensions, so that, for example, action is 

S = J dt dx 4 ~ £ C. 

After Hubbard- Stratonovich transformation, the Lagrangian density (2) takes the form 
C = ^ ^id t + ^pj* + & [p + p(x)]a 3 y 

- — 0*0 + ¥(T^<f)\ (3) 

c 

where is an auxiliary di-fermion field, \I/ is a two-component Nambu-Gor'kov field, ^ = 
(■0f,-0>J) T , o± = ± iffz), with o"i,2,3 being the Pauli matrices. In this calculation we 
will concentrate on the unitary regime, a = oo, which in the dimensional regularization 
corresponds to l/c = [1]. It is possible to use e expansion to describe system near 
unitarity [13], and the results of this paper may be extended to the near-unitary regime. 
This is left for future work. 

The density correlation function is defined by 

S(a;) = -i{n\T^(x)^(x) ^t(0M0)|n) logZ[p(x)]| p=0 , (4) 
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where \Q) is the ground state, 

Z\p(x)] = J D^D^ f D^D0*Exp i J C(ip,(j)-p(x)) 
is the generating functional and x = (t, x). 

III. LOW ENERGY BOSON PROPAGATOR OF THE SUPERFLUID PHASE 

Let us review some of the low-energy properties of the superfluid phase of the unitary 
Fermi gas which will be relevant for the density-density correlator calculation [10, 11]. 

For energies below the quasi particle mass, A ~ <f> , the only relevant degrees of freedom 
are the bosonic excitations of the superfluid condensate. The Lagrangian is given by 



/ 



W0) = -*TrLog| 2 |, (5) 

Po + g ip [x) % d t - ^ - /I 



where 0o is the uniform medium order parameter which to the NLO in e is given by 

0o = 2 hq [1 + 5(3 C - 1 + Log 2)] , C « 0.14424, /i = fi/e, (6) 
[1], while the coupling, 

m \ 2n J 

was chosen by Nishida and Son [1] in order to give bosons canonical kinetic term. 

Eq. (5) may be expanded in powers of and its derivatives as well as in powers of e. To 
NLO in e and to the second order in p , e p we have 

Left = l&V- 1 ®, (8) 

where <£> T = (0, 0*) and 

^n(p) = V 22 (-p)=Z(p --)+2» + - ^_-_ + _ 



e 2 ^-^2 / ^ n 3 t 3 
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^ W = u-i(p) = _£|» + £iE-ii-M + o £ J, £ |) , (9) 
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with p = (po,p), e p = p 2 /2m, Z = 1 + e(log 2 - 7)/2, where 7 « 0.57722, is the 
Euler-Mascheroni constant. To NLO in e the low-energy effective Lagrangian only con- 
tains quadratic terms in 0. The non-Gaussian terms are suppressed by higher powers 
of e. The situation is analogous to the chiral perturbation theory of the Color-Flavor- 
Locked phase of asymptotically dense quark matter where meson interaction terms are 
exponentially suppressed [14]. A generic term in the NLO Lagrangian is of the form 
&Cr S e pgep/A r+s_1 $, r + s > 1, c rs is 0(1); Po, e p < A are treated as 0(1) in e expansion. 

In the e expansion approach observables are calculated in the perturbation theory where 
fermions interact with bosonic excitations of the condensate with the coupling g oc \fe of (7). 
Detailed descriptions of Feynman rules at both zero and finite temperature are presented in 
[10, 13]. 

Superficially, to calculate density correlator to the NLO in e it should be sufficient to use 
the leading order (LO) boson propagator which is proportional to (p — e p /2) -1 . However, 
it does not have the superfluid mode pole at po = c s \p\ as |p| — > 0, where c s is the speed of 
sound, which we know should be present in the spectrum. It has been shown that to capture 
the right form of the superfluid mode excitation energy it is sufficient to use boson propagator 
with momentum independent 0(e) boson self-energy corrections resummed [10, 11]. Then 
the relevant NLO boson propagator is 

^>» = rf-4 + +*) + « -' = fft - (10) 

We will need sound velocity to NLO. From the NLO boson propagator (8) and (9) using 
(6) one gets 
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c„ = 



2m 



1 + \ 



(11) 



IV. THE DENSITY CORRELATION FUNCTION 

To NLO in e the density-density correlation function will consist of the contributions 
from Fig. (1): the one boson exchange diagram which produces 0(1/ e) contribution, and 
the 0(1) one fermion loop diagram. 

The dashed line in Fig. (1) is the dressed boson propagator discussed in the previous 
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FIG. 1: Feynman diagrams relevant for the calculation of the density correlation function to NLO 
in e. Thin solid lines depict propagating fermions in the superfluid background (lines with squares 
correspond to the anomalous propagators), dashed line is the dressed boson propagator, thick solid 
line is the source coupled to the particle density. The 0(l/e) diagram contributes to the imaginary 
part of the density-density correlator below quasiparticle threshold, lo < 2 A, where A ~ </>o is the 
quasiparticle mass. 

section (Eq. (10)). Thin solid line represents fermion propagator 



density. The coupling g is given in Eq. (7). For energies below twice the quasiparticle gap 
the only contribution to the imaginary part of density- density correlator will come from the 
one boson exchange diagram. 

The one boson exchange diagram and the one loop graph contributions are given by 






g 2 tr [G(p) a 3 G(p + k) <r_] V(k) u tr [G(q) a 3 G(q - k) o+] + 



+ (k 



<-> -k), 



(13) 



and 




(14) 
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respectively, and J k = J d d+1 k/ (2 7r) d+1 . The fermion loop integrals are calculated in d — A—e 
dimensions expanding to the second order in ho, e^, and a±, 03 were defined below Eq. (3). 
It was argued in [1] that fi was to be counted as 0(e). To the NLO /i ~ e is omitted 
everywhere in (13) except in the denominator of D(k) u . 

The result normalized to the NLO particle density, n, [1] is 

1 , , = 2e p (l-f)-^(6pg + eS))+0( e 2 , e^M) _ 
n [P ° ,ep) pl-e p (/i+ e f) + i8 + 0(ee p fi) 



£ 



+ 



+ 



where 



n 



My/* 4 



4^o V" 6/i ' 24 ' 120 
< Po < 2A, 5 = + , 



l + £ (6C-^ + 2Log2-^ 



M + O U si) 



(15) 



, 7s ^ 0.57722, 0.14424, (16) 



/x = /V £ ~ 00, and A//x — 2(1 — 0.345e) is the NLO quasiparticle mass [1], and the NLO 
expression for the uniform medium superfluid order parameter (6) has been used. The boson 
exchange contribution to (15) is not fully expanded in the powers of po, e p to preserve the 
pole structure. The dynamic structure factor, cr(u>,e p ), is, then 

' 2 ^( 1 -f)-6w(6^ + 4))' 



c(Po,e P ) = -bmS^Cp) = vr 



n 



2 Po 



x 



(17) 



where n is the NLO equilibrium density of Eq. (16). And the inelastic form factor (per unit 
volume), iS>me/(e p ), in this approximation is given by 

5e e 7 



<Si„ei{e p ) = — — / dulmS{uj,e p ) = 1 J— 



l- £ -- 

4 24/iQ 



(18) 



V. SUM RULE CHECKS 



Let us check the correlation function (15) against several identities for a density correlator. 
Using S(k) of (15) we find that 

1. the dispersion relation 



S(u = 0,e fc ) 



7T 



da; 



lmS(u,e k ) 



(19) 



is satisfied to the NLO in e and to 0(e 2 k ); 



2. the compressibility sum rule 



S(u = 0,e k = 0) 



n 



(20) 



mc 2 s 



where n is the equilibrium density given by (16) and c s is the speed of sound given 
by (11) is satisfied to the NLO in e. Note that the l.h.s. of (20) only contains the LO 
expression for c s . One also observes that the 1 — e/4 coefficient at e p in the numerator 
of the fraction in (15) is necessary to satisfy (20) to NLO; 

3. the energy weighted sum rule which for a unit volume of the medium is 



Let us emphasize that sum rule (21) is sensitive to the dynamics of quasi particles at 
uo > 2 A which we have not included in the calculation. We see that the sum rule 
is satisfied at LO in e. At NLO the e p coefficient is reduced to 1 — e/4 instead of 1 
prescribed by (21) which is due to the omission of the quasi particle contribution above 
the 2 A threshold. Also higher order terms suppressed by powers of e p //i ~ e p /A 
appear which is natural within the framework of the low energy effective theory. So, 
we see that including only the condensate excitation contribution to the NLO Im S(p) 
and omitting the rest of the (heavy) intermediate states we do not violate the energy 
weighted sum rule (21) too badly. 

These checks indicate that the density-density correlator given by Eq. (15) reproduces 
response properties of the medium at low energy reasonably well. 

VI. CONCLUSIONS 

We have computed density correlation function of unitary Fermi gas, S(u>, e p ), to NLO in 
the e expansion for energies below quasi particle threshold cu < 2A and quoted expression 
for low energy dynamic structure function and static structure factor. The simplest extrap- 
olation to d — 3 may be made by simply setting e — 1 in the NLO expressions (15,17,18). 




(21) 



for any e p . From S(cu, e p 



,) of (15) we get 




(22) 
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As is well known (see, e.g., [15]), dynamic structure function is proportional to the inelas- 
tic scattering cross-section of any external probe coupled to density, and may be measured 
experimentally. 

Based on the calculation one predicts a significant change in the dynamic structure func- 
tion at the quasiparticle threshold, u = 2 A, which may help experimental measurement of 
the quasi particle gap. 

Future work will include calculation of the dynamic structure function for energies above 
quasi particle threshold and of the spin response function relevant for the description of 
neutrino coupling to the dilute neutron matter (e.g., [16, 17]). 
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